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Abstract.  The  non-contact  acousto-thermal  signature  (NCATS)  is  a  nondestructive  evaluation  technique  with 
potential  to  detect  fatigue  in  materials  such  as  noisy  titanium  and  polymer  matrix  composites.  The  determination  of 
underlying  physical  mechanisms  and  properties  may  be  determined  by  parameter  estimation  via  nonlinear  regression. 

The  nonlinear  regression  analysis  formulation,  including  the  underlying  models,  is  discussed.  Several  models  and 
associated  data  analyses  are  given  along  with  the  assumptions  implicit  in  the  underlying  model.  The  results  of  these 
analyses  are  discussed. 


Background  and  Experimental  Procedure 

Non-Contact  Acousto-Thermal  Signature  (NCATS)  uses  the  conversion  of  acoustic  energy  to  heat  to  characterize 
evolving  damage  in  materials.  When  high  amplitude  acoustic  waves  encounter  a  material  boundary  they  both  reflect 
from  and  transmit  into  the  material.  The  energy  absorbed  in  the  material  is  subsequently  converted  to  heat  through 
the  diffusion  of  transverse  thermal  currents,  inter-crystalline  thermal  currents,  and  dislocation  motion  [1,  2].  In  these 
experiments,  high  amplitude  acoustic  waves  propagated  through  an  air  gap  and  interacted  with  a  Titanium  sample. 
The  experimental  setup  to  quantify  the  thermal  response  consisted  of  an  ultrasonic  horn  operating  at  20  kHz,  an  IR 
camera,  a  flat  specimen,  and  a  servo  hydraulic  machine  that  is  used  to  apply  fatigue  cycles  to  the  specimen.  The  20 
kHz  sonic  horn  was  placed  480  jam  from  the  specimen  as  shown  in  Figure  1(a)  and  was  pulsed  for  1000  ms.  The 
sonic  horn  does  not  contact  the  specimen  even  at  full  amplitude.  Temperature  rise  in  a  Titanium  specimen,  due  to  the 
conversion  of  sonic  energy  to  heat,  was  measured  by  an  Indigo-Merlin  Infrared  (IR)  camera  located  190  mm  from 
the  specimen  also  shown  in  Figures  1(a)  and  (b).  This  work  evaluated  experimental  data  from  2  trials  at  80%  sonic 
horn  power  and  3  trials  at  100%  power.  A  sophisticated  analysis  of  one  trial  at  100%  power,  including  curve  fitting,  is 
presented. 

A  custom  developed  system  and  software  were  used  to  control  the  sonic  horn  pulse  and  the  acquisition  of  ther¬ 
mal  response  data  from  the  IR  camera.  The  firing  of  the  sonic  horn  and  the  thermal  response  data  collection  were 
synchronized,  with  the  start  of  the  horn  pulse  equating  to  time  zero.  Thermal  response  data  was  collected  at  15  frames 
per  second  for  a  total  of  30  seconds,  and  data  was  simultaneously  collected  at  multiple  points  in  the  IR  camera  field 
of  view.  Temperature  data  used  in  this  analysis  begins  when  the  sonic  horn  pulse  ceases  at  1000  ms.  The  spatial 
resolution  of  the  IR  camera  for  this  experimental  configuration  was  determined  to  be  250  /am. 

The  Ti-6A1-4V  specimens  chosen  for  this  study  were  flat  rectangular  bars  having  a  width  of  50.8  mm,  a  thickness 
of  3.2  mm,  and  a  nominal  length  of  200  mm.  The  specimens  used  in  this  study  contained  a  geometric  feature,  a  6.4 
mm  diameter  open  hole.  The  Ti-6A1-4V  material  has  an  equiaxed  microstructure  with  an  average  grain  size  of  20  jam. 
To  increase  the  emissivity,  the  specimen  was  painted  black  in  the  region  of  interest  on  the  side  of  the  specimen  facing 
the  IR  camera. 


Data  Analysis 

The  initial  goal  of  the  data  analysis  was  to  identify,  detect  or  discern  the  physical  methods  occurring  in  the 
NCATS  experiment.  This  task  requires  sophisticated  data  analysis  techniques  including  hypothesis  testing  and  in¬ 
formation  theoretic  based  model  selection  (such  as  the  Akaike  Information  Criterion- AIC).  Both  hypothesis  testing 
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FIGURE  1.  (a)  NCATS  experimental  setup  with  sonic  horn,  specimen,  and  IR  camera  (b)  Sonic  Horn  and  IR  camera  spacing 
relative  to  the  specimen 


(see  [3,  4]  for  examples)  and  the  AIC  (an  example  is  given  in  [5]  and  a  thorough  introduction  may  be  found  in  [6]) 
require  that  parameter  estimation  be  carried  out.  In  order  to  carry  out  parameter  estimation,  a  model  (which  depends 
on  quantities  of  interest  and  the  independent  variables)  and  data  are  needed.  The  data  must  include  the  measurement 
of  the  signal  and  the  value  of  the  independent  variable  at  that  measurement.  Initial  data  analyses  were  carried  out  us¬ 
ing  single  location  data  (with  dependent  variable  time),  and  a  ID  thermal  model.  These  produced  anomalous  results. 
Thermal  models  which  capture  the  3D  geometry  of  the  specimens  were  considered.  In  order  to  analyze  the  ensemble 
data  of  the  36  cursor  locations  on  the  back  surface  over  all  the  times,  all  the  independent  variables  must  be  known.  The 
independent  variables  for  this  experiment  are  time  and  spatial  coordinates.  For  some  of  the  data  sets,  the  value  of  the 
spatial  coordinates  were  not  known.  Data  that  was  collected  previously  [2]  were  used  instead.  The  videos  were  used 
to  compare  with  the  specimen  to  determine  the  location  of  the  pixels  relative  to  each  other.  The  pixels  were  found  to 
be  0.25  mm  by  0.25  mm  which  confirmed  the  measurements  from  previous  work. 

In  the  classic  paper  by  Parker  et  al.  [7],  the  authors  derive  a  flash  method  to  determine  several  thermal  properties 
of  materials,  including  the  thermal  diffusivity.  The  derivation  of  this  method  is  discussed  below.  Throughout  the 
discussion  the  method  and  estimate  are  referred  to  as  the  “Parker  method”  and  “Parker  estimate,”  respectively.  The 
Parker  estimate  of  the  thermal  diffusivity  ap  is  obtained  by  making  several  assumptions  about  the  sample.  Using  these 
assumptions,  a  series  solution  may  be  derived  and  then  truncated  to  obtain  the  Parker  result.  The  assumptions  implicit 
in  using  this  method  are  equivalent  to  assuming  the  temperature  u  is  the  solution  to  the  heat  equation  which  is  given 
by  the  partial  differential  equation 

ut  -  auzz  =  0  for  0  <  z  <  L.  (1) 

The  thermal  diffusivity  is  given  by  a  mm2s-1  and  the  length  that  the  heat  passes  through  is  given  by  L  mm.  Parker’s 
method  further  assumes  no  loss  on  the  boundaries  (perfectly  insulated  boundaries)  which  corresponds  to  Neumann 
boundary  conditions 


-auz\z=0  =  0  (2) 

auz\z=L  =  0. 


In  Parker’s  method  for  flash  heating,  rather  than  representing  the  flash  by  a  source  on  the  boundary,  the  flash  is 
assumed  to  be  uniformly  absorbed  in  a  small  region  near  the  heated  boundary  (the  z  =  0  boundary  here)  yielding  the 
Heaviside  initial  temperature 


u(  0,  z) 


tl  =  — 

^  PCpg 

0 


for  0  <  z  <  g 
for  g  <  z  <  L 


(3) 


where  Q  (J  cm  2)  is  the  radiant  energy  of  the  flash,  g  (mm)  is  the  depth  that  it  is  uniformly  absorbed,  p  (g  mm  3)  is 
the  density  and  cp  (J  kg_1K_1 )  is  the  heat  capacity.  A  well  established  result  from  Parker  is  that  given  the  temperature 
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Parker  fit 


(a) 


(b) 


FIGURE  2.  (a)  Parker  estimates  of  the  thermal  diffusivity  at  each  location,  for  each  trial,  (b)  An  example  fit. 


corresponds  to  a  solution  of  equations  (l)-(3),  the  thermal  diffusivity  may  be  estimated  as 

&P  =  1.38L2/A/2,  (4) 

where  A/2  is  the  time  that  the  back  surface  reaches  1/2  of  the  maximum  temperature.  The  estimates  for  each  location 
on  the  sample  are  given  in  Figure  2(a).  The  handbook  value  for  the  thermal  diffusivity  of  Ti-6A1-4V  is  2.75  mm2s-1 
[8].  The  Parker  estimates  obtained  by  applying  Equation  4  to  the  data  from  each  of  the  36  locations  from  each  of  the 
five  trials  are  depicted  in  Figure  2(a).  It  is  clear  in  this  figure  that  the  Parker  estimates  are  about  half  the  handbook 
value.  A  further  problem  with  the  Parker  estimates  is  that  the  implicit  assumptions  are  clearly  not  valid.  One  clear 
violation  of  these  assumptions  is  loss  on  the  boundary.  Cooling  may  be  observed  after  two  seconds  in  the  data  from 
one  location  during  one  trial  (see  Figure  2(b)).  This  model  assumption  can  be  relaxed  by  replacing  the  insulated 
boundary  condition  in  Equation  3  with  boundary  conditions  with  losses 


-auz\z=0  =  -A0u(t,  0)  (5) 

auz\z=L  =  -Aiuit.L)  (6) 

where  To  and  Al  (mm  s-1)  are  the  losses  on  the  right  and  left  boundary  respectively.  This  model  was  approximated 
numerically  by  the  finite  element  method  [9].  The  left  side  temperature  TL-  and  the  absorption  depth  g  in  Parker 
initial  condition  in  (3),  the  thermal  diffusivity  in  the  heat  equation  (1)  and  the  loss  terms  in  the  boundary  conditions 
(6)  were  estimated  using  nonlinear  least  squares.  The  estimated  thermal  diffusivity  was  0.225  mm2s_1  which  is  an 
order  of  magnitude  less  than  the  handbook  value  of  2.75  mm2s_1.  The  resulting  model  values  also  do  not  fit  that  data 
well  as  can  be  seen  in  Figure  2(b). 

The  poor  fit  and  lack  of  correspondence  of  the  estimated  thermal  diffusivity  with  the  handbook  value  indicate  that 
there  may  be  other  invalid  assumptions  that  are  effecting  the  results.  The  thermal  image  of  the  specimen  12  ms  after 

the  horn  is  turned  off  in  Figure  3(a)  suggests  that  the  heating  is  nonuniform.  There  does  appear  to  be  radial  symmetry 

which  leads  us  to  consider  a  radially  symmetric  cylindrical  domain.  It  is  further  assumed  that  the  edges  are  far  enough 
away  to  be  considered  infinite.  The  corresponding  partial  differential  equation  is  given  by 


ut  -  a 


=  0 


^rlr=0  —  0 


lim  u  =  0 

r— »oo 

-auz\z=0  =  -Aou(t,  r,  0) 


(7) 
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Max  Temperature 


FIGURE  3.  (a)  An  example  thermal  image  from  the  experiment  (b)  The  maximum  temperature  sensitivity  to  amplitude 


auz\z=L  =  -Aiu(t,r,L) 
u( 0,  r,  z)  =  uo(r,  z ) 

where  r  is  the  radius  and  z  is  the  depth  in  the  specimen  (in  mm).  For  the  cylindrical  parameter  estimation,  the  losses 
on  the  boundaries  were  assumed  to  be  equal,  so  To  =  Ai  =  A.  The  partial  differential  equation  with  boundary  and 
initial  conditions  (7)  was  solved  with  finite  differences  and  the  finite  volume  method.  The  results  presented  here  were 
computed  using  the  finite  volume  method. 

The  logistic  function  is  a  smooth  approximation  of  the  Heaviside  function.  For  instance,  the  Parker  initial  condi¬ 
tion  (3)  may  be  smoothly  approximated  by 

^  =  1  +  el5iz-g)  ’  ® 

where  the  constant  15  was  chosen  arbitrarily.  The  radially  symmetric,  with  discontinuity  at  (r,  z)  =  (ro,  Zo)  rather  than 


FIGURE  4.  (a)  Estimated  initial  temperature  distribution  using  (9)  (b)  Estimated  fit 
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z  =  g  as  in  (8),  version  of  this  is  given  by 


wo  (r,z) 


_ Tl _ 

(1  +  e15^~z^)(l  +  e15^)) 


+  Tr , 


(9) 


the  right  side  temperature  TR  is  fixed  to  be  the  value  of  the  observed  temperature  at  t  -  Is.  This  time  value  was  chosen 
because  it  is  when  the  horn  was  turned  off.  The  subsequent  data  that  is  presented  here  begins  at  t  >  1.  The  above 
initial  condition  gives  a  temperature  profile  that  is  approximately  Tl  for  depths  less  than  zo  (z  <  zo)  and  radii  less  than 
r0  (r  <  ro),  and  TR  everywhere  else.  Using  this  initial  condition,  the  estimated  thermal  diffusivity  (cr)  and  loss  were 
approximately  0.22.  The  estimated  thermal  diffusivity  was  is  an  order  of  magnitude  less  than  the  handbook  value.  The 
estimated  initial  condition  (Tl,  Zo,  and  r0  were  estimated)  is  given  in  Figure  4(a).  The  curve  fit  using  this  model  and 
initial  condition  were  much  better  than  the  fit  using  the  ID  models.  The  curve  fit  is  given  in  Figure  4(b). 


Initial  Temperature  Distribution 


Model  Fit 


(a) 


(b) 


FIGURE  5.  (a)  Estimated  initial  temperature  distribution  using  (10)  (b)  Resulting  fit 


The  initial  condition  (9)  when  used  in  data  analysis  yielded  anomalous  thermal  diffusivity  estimates.  Splines 
were  then  used  to  estimate  the  z-dependence  of  the  initial  temperature  distribution.  The  initial  condition  was  given  by 


uo(r,z)  = 


TLf(z) 

(1  +  ^15(r-r0))  R 


(10) 


where  /(z)  is  a  cubic  polynomial  with  nodes  at  z  =  0  and  z  =  L.  The  value  at  z  =  L  was  fixed  to  be  zero,  the 
temperature  at  z  =  0,  and  derivatives  with  respect  to  z  at  the  nodes  were  estimated.  The  estimated  initial  temperature 
distribution  and  model  fit  are  depicted  in  Figures  5(a)  and  (b),  respectively.  The  results  were  similar  to  fit  using  the 
logistic  dependence  in  z.  The  similarity  in  the  estimated  curves  is  not  surprising.  Estimating  the  initial  temperature 
distribution  in  the  heat  equation  is  ill-posed  in  the  sense  that  many  different  initial  temperature  profiles  will  yield 
the  same  or  similar  results.  This  can  be  seen  by  comparing  the  temperature  profiles  in  Figures  4  and  5(a)  and  then 
comparing  the  resultant  models  in  Figures  4  and  5(b).  The  estimated  thermal  diffusivity  was  also  an  order  of  magnitude 
lower  than  the  handbook  value.  When  another  node  was  added  to  the  spline,  the  results  were  similar. 

Many  other  initial  temperature  distributions  and  boundary  conditions  were  used.  A  trigonometric  source  was 
used,  with  a  poor  fit  and  thermal  diffusivity  within  the  same  range  as  the  other  thermal  diffusivity  estimates.  Initial 
temperature  distributions  that  are  not  presented  here  were  also  used,  including  products  of  trigonometric  and  expo¬ 
nential  and  products  of  splines  and  exponentials.  These  fits  were  also  poor  with  thermal  diffusivity  estimates  that  were 
an  order  of  magnitude  lower  than  the  handbook  value  of  the  thermal  diffusivity. 


Conclusions  and  Future  Work 

The  work  here  illustrates  the  difficulties  with  data  analysis.  The  results  were  anomalous.  The  observed  thermal 
diffusivity  estimates  were  not  within  the  range  of  the  handbook  values  though  this  behavior  is  consistent  with  these 
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experiments  carried  out  on  other  Ti-6A1-4V  samples  with  different  geometries.  The  sample  will  be  interrogated  with 
other  methods  to  determine  its  thermal  diffusivity.  The  behavior  of  the  horn  and  possible  other  physics  will  also  be 
investigated.  Further  work,  once  the  estimation  for  a  single  location  is  carried  out,  will  be  done  to  incorporate  all  of 
the  observed  locations. 
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